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Abstract 

This paper evaluates the results of a microprocessor based fault locator employed on a laboratory power system model with 
those obtained within EMTP/ATP simulations. The Discrete Fourier Transform concept is used to extract the fundamental 
frequency components of the voltage and current waveforms under fault and the reactance based method to locate the fault. 
The results suggest a close correlation between the experimental and simulation studies for all fault types. 
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Introduction 


Distribution line faults have been a primary concern to networks. Faults affect the reliability of the service to the 
customer, the stability of the power system and the quality of the power delivered. Distribution line faults can even 
cause complete supply outages. These faults must be detected quickly for fast deactivation of the faulty section of 
line, in order to reroute the power to the customer to meet target restoration times. 


Distribution line fault locators can be classified under two categories. The first is the impedance based locators, 
which analyse the power frequency (50 Hz or 60 Hz). The second is the travelling wave based locators which 
analyse the travel time of high frequency transients. Transient based fault location methods require higher 
sampling rates and more responsive voltage or current transducers. They also require multiple fault transient 
detectors placed throughout the system, which are clocked using GPS technology and data sent to a central 
location via a communication link. The use of GPS can make such a solution very expensive. Transient based fault 
location techniques for radial, meshed and distribution lines with tapped load power system configurations are 
described by a number of researchers (Nouri et al 2014, Magnago et al 1998 and Nouri et al 2001), respectively. In 
impedance based methods, the ratio of extracted fundamental frequency of the voltage and current waveforms 
yields the impedance seen from the generator up to the point of fault. A useful summary of the available 
impedance calculations can be found in (Zimmerman 2005). The reactance equation that will be used to obtain the 
fault location (m) in per unit is expressed as: 



a) 


where IM defined as storing of the imaginary part of the voltage ( v v ) and the current ( i s ) ratio at the sending end 
of the line, and Xitl is the positive sequence reactance of the distribution line. 


Discrete Fourier Transform 

Mathematical Background 

In power systems the voltage and current waveforms are essentially shifted sine or cosine waves of fundamental 
frequency. These signals are effectively periodic when the system is running under steady state conditions. Under 
fault conditions, the fundamental frequency of the voltage and current amplitude can be extracted from the non 
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periodic signals using the Fourier Transform. The basis functions chosen for the Fourier Transform is the complex 
exponential g(/) = e~ J(Dt . The Fourier Transform is defined as follows: 

00 

F(a>)= \f(t)e ia dt ( 2 ) 

—oo 

On closer examination of (2), it can be seen that Fourier Transform integration takes place over all of time, i.e. from 
-oo to oo . Therefore, the Fourier Transform should be reserved for periodic signals only. A power system 
contains many devices such as transformers, loads and switching circuits. These all create harmonic distortion 
under steady state conditions. The Fourier Transform is useful to analyse the harmonic content of the signal. 
However, when the system is running under a fault condition, non periodic signals exist. Consequently, the 
Fourier Transform should not be used to analyse fault conditions. An alternative form of Fourier's analysis 
techniques should be used, known as the Fourier Series. Unlike the Fourier Transform, the Fourier Series analyses 
the signal over a finite period of time. The basic assumption used in Fourier Series is that the waveform that results 
from the fault condition is assumed to be periodic within the interval from t 0 to t 0 + T , where T is the period of 
the fundamental component. The following equations are derived from (Wood 1984, PhadKe et al 1998, Johns et al 
1995, and Croft et al 2001). 

The voltage or current waveform s(t ) can be represented using the Fourier Series. 

a () 00 00 

s(t ) = — + X a n cos nco 0 t + X b n sin nco 0 t (3) 

2 n= 1 n = 1 


where the Fourier coefficients are: 

2 to+T 

a n = — J sj/ 1 ) cos nco 0 t dt 

T to 

And 


(4) 


2 

b n - — J s(t) sin neo 0 t dt 


(5) 


where n = 1, 2 , 3 ... , for the fundamental, second, and third harmonic ... . It is apparent that the fundamental 
frequency can be selected by setting n- 1 . Thus, the Fourier coefficients for the fundamental frequency can be 
determined: 

2 to+T 

a x -— \ s(t ) cos co 0 t dt ( 6 ) 

T to 

And 


by 


2 _ 

T 


L+T 


J s(t) sin co 0 t dt 


(7) 


In discrete analysis the time of the jth sample is tj - jAt and the period of the waveform is T = N At , where N is the 

number of samples per cycle of 50Hz. Now ( 6 ) and (7) can be determined in discrete time: 

2 n- 1 2 N ~ l 2 N ~ l 

<h = ~777~ £ s(t ) cos co 0 t ■ = — — X s(j At) coso) JAt = — X s(j) cos 0 ) JAt ( 8 ) 

NAt j = 0 NAt y=o N 7=0 


And 


2 N - 1 2 N-l 2 N-l 

by = T7— Z s(tj ) sin O) 0 tj = —— X s(jAt) sin coJAt = — X s(j) sm coJAt (9) 

NAt 7=0 NAt 7=0 N 7=0 

The amplitude of the real and imaginary parts of the fundamental frequency signal are represented by the Fourier 
coefficients a x and b x respectively. It is necessary to determine the amplitude of the fundamental frequency in 
order to apply an impedance equation to locate the fault. 
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s 



( 10 ) 


Also it is required to find the phase angle between the real and imaginary components of the fundamental 
frequency. 


Zs = tan 


r 7 'n 


(ii) 


Practical Application of DPT 

A single line diagram of the laboratory power system model is shown in Fig. 1, with a fault locator system 
connected to bus-A. It can be assumed that the distribution line has a length of 1 pu. 

A system must be designed to gather the phase voltages and currents and upload them to the CPU to analyse the 
fault conditions. This requires a sensor circuit to measure the voltages and currents, interfaces between the 
analogue and digital domains, and an external microprocessor to control the data gathering process. The chosen 
external microprocessor is the Siemens C164CI on the PHYTEC kitCON-164 evaluation board. A block diagram of 
the required system is shown in Fig. 2. The ADC and DAC are 12 bit devices. The data sent between the C164CI 
microcontroller and the CPU is via an RS-232 bus and the method for the C164CI communicating with the ADC 
and DAC is a serial peripheral interface (SPI) bus. 

A photograph of the laboratory power system model and microprocessor based fault locator is shown in Fig. 3. 



Fault Locator 
System 


FIG. 1 SINGLE LINE DIAGRAM OF THE LABORATORY POWER SYSTEM MODEL 

A 



FIG. 2 A BLOCK DIAGRAM OF THE DATA ACQUISITION SYSTEM 



FIG. 3 PHOTOGRAPH OF LABORATORY POWER SYSTEM MODEL AND MICROPROCESSOR FAULT LOCATOR 
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FIG. 4 BLOCK DIAGRAM OF SOFTWARE FUNCTIONALITY 


Fault Location Software 

The fault location software is a function based program as shown in Fig. 4. There are three main functions which 
determine the fault location: - 

i. The function getsamples() uses timer 2 interrupt service subroutine to generate a sampling period of 625 ps. 
At each sampling interval the function collects the instantaneous value of the line voltages and currents. It repeats 
this until the required number of samples is gathered. 

ii. This data is passed on to the function DFT(), which extracts the amplitude and phase angle of the 
fundamental frequency. 

iii. The analyse() function uses a fault current threshold determined during program initialisation to decide 
which phases are faulty. A variation of reactance equation is then applied to locate the fault depending on the 
number of faulty phases. The reactance equation gives a fault location in mpu. 

EMTP-ATP Simulation 

A line-to-ground on phase A is simulated using an EMTP/ATP simulation model shown in Fig. 5. The fault occurs 
in the middle of the distribution line with a fault location of 0.5 pu. The EMTP/ATP model parameters are based on 
the laboratory power system model shown in Fig. 1. 



FIG. 5 ATP SIMULATION MODEL OF A FAULTY 415 V DISTRIBUTION LINE 


Results & Discussion 

EMTP-ATP Simulation 

A fault is set to occur at t = 0.06 seconds and lasts for the remaining duration of the signal (5 cycles). The profiles of 
the ATP voltage and current waveforms can be seen for phases A, B and C in Fig. 6 and Fig. 7, respectively. 

The DFT is used to extract the amplitude and phase angle of fundamental frequency line voltages and currents 
from Fig. 6 and Fig. 7. All transmission line signals are sampled at a rate of 32 samples per 50 Hz cycle, which 
produces a sampling period of T s = 625 jus . The DFT fundamental frequency extraction of phase A for voltage and 
current can be seen in Fig. 8 and figure Fig. 9, respectively. By examining these waveforms it can be seen that there 
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are 64 pre-fault samples (-64 to 0) and 160 samples after the time of fault occurrence. During the pre-fault samples, 
the fundamental frequency voltage and current amplitude is 1 pu. The application of a fault on the line causes the 
current at the sending end to increase to approximately 11 pu. Also, there is a decaying exponential DC offset of 
short duration. It should be noted that the DFT ignores this DC offset and just tracks the amplitude of the 
fundamental frequency, as required. The fundamental frequency voltage amplitude does not change by a 
significant amount. 



FIG. 6 ATP PHASE VOLTAGES WAVEFORMS FIG. 7 ATP PHASE CURRENTS WAVEFORMS 



FIG. 8 EMTP-ATP FUNDAMENTAL FREQUENCY VOLTAGE FIG. 9 EMTP-ATP FUNDAMENTAL FREQUENCY CURRENT 
EXTRACTION USING DFT EXTRACTION USING DFT 


These results show that the DFT can extract the fundamental frequency in approximately 32 samples, i.e. a 
complete cycle of 50 Hz. The resultant fault location in pu can be found by using the reactance equation (discussed 
earlier) as shown in Fig. 10. 

This waveform correctly shows that the fault location converges to 0.5 pu after the application of the fault. 

The resultant fault location in per unit can be found using the reactance equation as show in Fig. 13. 



FIG. 10 EMTP-ATP FAULT LOCATION WAVEFORM FIG. 11 MICROPROCESSOR BASED FUNDAMENTAL 

FREQUENCY VOLTAGE EXTRACTION USING DFT 
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FIG. 12 MICROPROCESSOR BASED FUNDAMENTAL FIG. 13 MICROPROCESSOR FAULT LOCATION WAVEFORM 

FREQUENCY CURRENT EXTRACTION USING DFT 

Verification of Results Using the Microprocessor Based Fault Locator 

Results produced by EMTP/ATP simulation are verified using the microprocessor based fault locator on the 
laboratory power system model. The DFT is used to extract the fundamental frequency, amplitude and phase angle 
for all three line voltages and currents. Similarly to the EMTP/ATP simulation a sampling condition of 32 samples 
per 50 Hz cycle is used. The DFT fundamental frequency extraction of the voltage and current waveforms can be 
seen in Fig. 11 and Fig. 12, respectively. 

Conclusions 

The results produced in the EMTP/ATP environment are compared with experimental results produced by 
microprocessor based fault locator on a laboratory power system model. The DFT extracts the fundamental 
frequency in less than a complete cycle of 50Hz and when combined with the reactance equation produces very 
accurate fault location. The DFT effectively ignores DC offset especially when present in heavily inductive lines. 

A collapse in line voltage or an increase in line current is a useful indicator of a faulty phase. The fault locator 
allows quick detection for fast deactivation of the faulty line. This allows the power to be rerouted to the customer 
to meet target restoration times. Typically, distribution line faults have a fault impedance of unknown quantity, 
causing error in fault location. The reactance method significantly reduces this error to a tolerable amount. 
Generally, fault impedance is resistive in nature. The reactance equation ignores all resistances in the system and 
locates the fault purely based on reactive impedance. 

The DFT requires a large amount of processing power due to the need for it to analyse a complete 50 Hz cycle at 
every sample period. It has the major advantage over other methods of extracting the fundamental frequency that 
no prior distribution line information is required apart from the reactive part of the line impedance. 
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